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ABSTRACT 

We present the results of a high angular resolution study of the BL Lac object Markarian 
501 in the radio band. We consider data taken at 14 different epochs, ranging between 1.6 GHz 
and 22 GHz in frequency, and including new Space VLBI observations obtained on 2001 March 
5 and 6 at 1.6 and 5 GHz. We study the kinematics of the parsec-scale jet and estimate its bulk 
yelocity and orientation with respect to the hne of sight. Limb brightened structure in the jet 
is clearly visible in our data and we discuss its possible origin in terms of velocity gradicnts in 
the jet. Quasi-simultaneous multi-wavelength observations allow us to map the spectral index 
distribution and to compare it to the jet morphology. Pinally, we estimate the physical parameters 
of the parsec-scale jet. 
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1. INTRODUCTION 

BL Lac objects are one of the several flavors 
of radio loud active galactic nuclei. Along with 
peculiar properties at other wavelengths (lack of 
strong emission lines, high levels of variability, op- 
tical polarization up to 3%, detection at X-ray and 
7-ray energies), they also show quite extreme be- 
havior in the radio: the parsec-scale structure is 
dominated by a compact, Aat-spectrum core, from 
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which a one-sided jet emerges, showing knots of 
enhanced brightness. When monitored over long 
intervals, these components are often found to be 
in motion, sometimes with an apparently superlu- 
minal velocity (see, e.g., Homan et al. 2001). 

The currently accepted explanation for these 
properties is found in the frame of unified models 
(see e.g. Urry & Padovani 1995). These models are 
based on a dust enshrouded super-massive black 
hole accreting matter in a disk, from which two 
symmetric collimated jets of relativistic plasma 
are ejected in opposite directions. When the angle 
between the jet axis and the line of sight is small, 
the resulting Doppler boosting accounts for many 
of the above-mentioned properties. 

A source such as Markarian 501 (B1652+39) is 
an ideal target to test the assumptions involved in 
unihed models. At its low redshitt {z = 0.034) ^*', 
Imas — 0.72 pc; therctorc, thc milliarcsecond res- 
olution achievable with VLBI tcchniques makes it 
possible to investigate in great detail the regions 
near the central core. Furthermore, Mkn 501 is 
one of the few sources with a clear limb-brightened 
structure (Giovannini 1999; Aaron 1999). Pinally, 
this source is very well studied at other frequencies 
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and it is known for its X— and 7— ray (Quinn et 
al. 1996; Bradbury et al. 1997) activity, which also 
sets constraints on the parameters describing the 
physics of the inner jet (sec Tavccchio ct al. 2001; 
Katarzyhski ct al. 2001, and relerences therein). 

In this paper we present a large multi-epoch, 
multi-frequency data-set of VLBI images; new ob- 
servations as well as previously pubhshed data 
are used to study the physical properties of this 
source. In particular, we will discuss the pc-scale 
jet morphology, velocity and orientation. Addi- 
tionally, our data-set allows us to compare images 
obtained at different frequencies with similar reso- 
lutions and thus to perform spectral index^^ stud- 
ies for comparison with the total intensity images. 

In § 2 we describe the previous and new obser- 
vations used in this paper, along with thc data 
reduction methods and procedures for component 
titting and spectral index mapping. Results are 
presented in § 3 and discussed in § 4. We present 
our conclusions in § 5. 

2. OBSERYATIONS 

Tablc 1 summarizcs all obscrvations considcrcd 
in this paper, listing epoch, frequency, total ob- 
serving time and array used. The notes provide 
references for dctails. Thc prcscnt data-sct is vcry 
similar to the one used by Edwards & Piner (2002, 
hereafter EP02) to discuss the proper motion in 
Mkn 501: wc did not use the 1995 Oct 17 cpoch 
(due to its low sensitivity), or the 1996 Jul 10 and 
1998 Oct 30 epochs (as better data from other 
ncarby cpochs was available). On thc othcr hand, 
we add 3 more epochs (see note 4 to Table 1) at 15 
and 22 GHz and we consider two new Space VLBI 
(VSOP) obscrvations obtained on 2001 March 5 
and 6 (see Table 2). 

2.1. Data Reduction 

2.1.1. Ground Observations 

Imagcs from somc of thc prcvious obscrvations 
have been published or made available to the sci- 
entiiic community on open web sites. However, for 
thc sakc of complctcncss and homogcncity, wc re- 
quested and obtained the calibrated (it, t^^-data for 
all observations presented here. After the initial 

We will define the spectral index a such that S^i/) oc 



reduction (scc notcs to Tablc 1 for references), we 
imported the calibrated data into the NRAO As- 
tronomical Image Processing System (AIPS) and 
inspectcd thcm caretuUy. Thcn, wc produccd ncw 
images; one or more self-calibration cycles proved 
to be useful in some cases, providing a better sig- 
nal to noisc ratio. Wc adoptcd thc appropriatc 
{u,v)-ra,nge, weights, cell-size and angular reso- 
lution to map the spectral index and to look for 
propcr motion of thc jct substnictures; we present 
parameters of hnal images in Table 3. 

2.1.2. Space VLBI Obseruations 

We obtaincd Space VLBI images at three dif- 
ferent epochs: 1997.8, 1998.4, and 2001.3. VSOP 
(VLBI Space Observatory Programme) observa- 
tions combinc data from ground arrays, such as 
the VLBA, and the 8m orbiting antenna on board 
the satellite HALCA. For the latter two epochs ob- 
servations at 1.6 and 4.8 GHz wcrc pertormed on 
successive days to minimize flux density variabil- 
ity problems in spectral index mapping. The first 
Space VLBI obscrvation was conducted at 1.6 GHz 
only. Eurther details on these observations are 
given in Table 2. 

Thc standard VSOP obscrving modc has one 
polarization (LCP) and two lFs of 16 MHz each 
with two-bit sampled data. The lowest IF lailed 
in the 1.6 GHz obscrvation in 1997. In 2001, the 
"a priori" calibration was affected by lack of sys- 
tem temperature measurement for Robledo and 
Goldstonc. A mini-cal proccdurc was donc bc- 
fore starting the experiment in Robledo, provid- 
ing a value of 20 K that was assumed as constant. 
For HALCA, wc uscd thc nominal systcm tcm- 
peratures of 75 K and 77 K for the two IFs at 
1.6 GHz, and 88 K and 92 K for those at 4.8 GHz 
(Hirabayashi et al. 1998, 2000). 

AU the data were correlated in Socorro and 
then imported and reduced in AIPS. Global fringe- 
titting returnc;(l good solutions for all the observa- 
tions in which large ground telescopes were used, 
although no solutions were found for the last ob- 
servation (2001 Mar 6), unless a valuc of SNR as 
low as 4.5 was set in the flt. In spite of such a low 
signal to noise ratio, we found these solutions rea- 
sonable: delays and rates found in this way do not 
present any discontinuity or suspicious behavior. 
Thus, we applied these solutions and proceeded 
with the whole array. 
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Wc thcn pcrformcd sclf calibration using clcan- 
component models, first in phase and, afterward, 
in both phase and amphtude. In the latter, we re- 
pcatcd thc task scvcral timcs, incrcasing thc base- 
line range. This revealed to be of great impor- 
tance for those observations where gain informa- 
tion for somc antcnnas was missing. For cxamplc, 
we could self-calibrate the shortest space-ground 
baselines using 'uniform weight' ground-only im- 
agcs; in fact, wc havc similar baselincs both in the 
ground array and to HALCA, as provided by its 
highly elliptical orbit (see Fig. 1, 2). Figures 3 
- 9 prcscnt all thc images, in order of decreasing 
angular resolution. 

2.2. Component Pitting 

Identihcation of componcnts is usually thc most 
ditScult problem affecting proper motion studies. 
Long time intervals between epochs can lead to 
misidcntification of componcnts; morcovcr, thc use 
of different resolutions can cause confusion, for 
example by blending components with different 
propertics and vclocitics; finally, rcgistcring im- 
ages at different observing frequencies is also nec- 
essary: the apparent core shift due to the different 
spectral index of the core and jet can be signifi- 
cant. 

In the present work, the time coverage is very 
dense (14 epochs ovcr 4 ycars, plus thc hnal VSOP 
observation) and there are no large time gaps. To 
prevent angular resolution problems, we produced 
maps at thc same resolution at all cpochs. Evcn 
if some observations allowed a better resolution, 
we convolved all the images obtained at a fre- 
qucncy 5 GHz < ;^ < 15 GHz with a circular rcstor- 
ing beam of 1.2 mas FWHM. However, higher res- 
olution maps have also been made, in order to pro- 
vide morc inlormation whcn nccdcd. In particular 
we obtained images at a resolution of 0.6 x 0.9 
mas (RA x Dec) in all observations with good 
(u, u)-covcragc on long basclincs (scc Tablc 3). A 
comparison of fit components in convolved and full 
resolution images revealed good agreement. 

The fit to the components of thc jct was done 
within AIPS and the Difmap package^^. In AIPS 
we used the task JMFIT, which allows one to fit 
Gaussian components to part of an image. Up 
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to four componcnts can bc used at a timc and 
the task, starting from given initial conditions, 
solves in an iterative way for position, flux den- 
sity and cxtcnsion of cach componcnt. Sincc more 
than four components are needed to properly fit 
the images, we had to run the task several times 
on diffcrcnt rcgions. Thc rcgions wcrc madc to 
overlap with each other in order to obtain results 
compatible with previously published models. We 
also rcpcatcd thc Ĕt looking for "stablc" solutions, 
i.e., changing our initial guesses or even leaving 
them blank. Eventually, very good models were 
obtained; thc diffcrcncc in flux dcnsity bctwccn 
our fits and thc images is always smaller than 5%. 

Using Difmap, we also carried out model-fitting 
to the visibilities rather than the image. Similarly 
to the fits made in AIPS, four to six jet compo- 
nents were required in addition to the core to pro- 
vide a good representation of the data. 

2.3. Spectral Index Mapping 

When producing spcctral index images, there 
arc thrcc main cffc!cts that havc to bc carc!fully 
taken into account. First, because of fiux density 
variability and possible proper motion of compo- 
ncnts, thc obscrvations must bc takcn at vcry close 
epochs. Next, the angular resolution must be the 
same at both frequencies and the (u, w)-coverage 
as similar as possiblc; this guarantccs thc same 
sensitivity to extended structures in the high and 
low frequency images. Pinally, as absolute posi- 
tion is not availablc, wc nccd to properly align the 
images before combining them. 

To meet the first two requirements, we pro- 
cluced spectral index images using data separated 
by less than 24hrs. We considcrcd 15 and 22 GHz 
data obtained on the same day (1997 Aug 15) and 
Space VLBI obscrvations at 1.6 and 5GHz made 
on successive days (1998 Apr 7 and 8 — sec Ed- 
wards et al. 2000, for a preliminary analysis). We 
note that 1998 Spacc VLBI observations have bet- 
ter (w, t;)-coverage than those in 2001 (see Fig. 1, 
2). To obtain matching (u, w)-covcrages, we cut 
the shortest baselines from the low frequency data 
and the longest baselines in the high frequency ob- 
servations. Images were then obtained using uni- 
form weights, and the same cell-size and restoring 
beam. 

The problem in registering the images is that 
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thc pcak of thc imagc, coincident with the core, 
can correspond to diSPerent positions at different 
£requencies if the core is self-absorbed. So, we 
ahgncd thc iniagcs using thc positions of compo- 
nents along the jet in a region with steep spectral 
index and no self-absorption effects. We found 
that at 15 and 22 GHz no shift was ncccssary, while 
a shift of ~0.2 mas was necessary in both RA and 
Dec between the 5 and 1.6 GHz images. See Fig- 
ures 10 - 12 for the iinal spectral index maps. 

3. RESULTS 

3.1. Source Morphology 

Thc largc scale radio structurc has bccn imagcd 
with the VLA (Ulvcstad et al. 1983; Van Brcugel 
& Schilizzi 1986; Cassaro ct al. 1999), showing a 
corc dominated sourcc with two-sidcd diffuse emis- 
sion oriented at PA ~ 45°; this indicates that the 
jets are non relativistic at a (projected) distance 
of ~ lOkpc from the core. The pc-scale structure 
has also been investigatcd, using VLBI tcchniques 
(see e.g. Conway & Wrobel 1995), with resolution 
of ~ 10 mas. However, the coUection of multi- 
epoch, multi-frequency, high resolution and high 
sensitivity images will improve our knowledge of 
the source morphology and evolution. The present 
VLBI observations show a strong core and a one- 
sided jet. Thc jet exhibits multiple sharp bends 
before undergoing a last turn, followed by rapid 
expansion. We can distinguish three different re- 
gions: 

• a first region, extending ~ lOmas from the 
core, whcre a high brightness jet structure 
is present (Fig. 3, 4, and 5). The jet PA 
is not constant in this region, being 150°- 
160° near the core (< 1 mas) and moving to 
^ 145° (£rom 1 to ~ 4 mas from the core) to 
become ~ 90° from 4 to 10 mas. 

Thc 22 GHz high rcsolution imagc (Fig. 3) 
shows that the jet is resolved and limb- 
brightened starting at ~ 1 mas from the core. 
This is conhrmcd by thc high scnsitivity 
5GHz Space VLBI image (Fig. 5); although 
this image does not have enough angular res- 
olution to resolve the inncr jct ncar to thc 
core, it shows an extended low brightness 
emission, contirming the transverse exten- 
sion of the jet in the inner 2mas. In the 



rcgion from 2 to lOmas from thc corc, thc 
limb brightened jet structure is clearly visi- 
ble. 

The brightness proAle along the main axis 

of the jet is not unitorm; high brightncss rc- 
gions are followed by dimmer regions. In 
particular at 5-6 mas from the core a dccp 
minimum is prcscnt just betore an extcnded 
bright spot with a pcak 8-9 mas from the 
core. However, we do not refer to thcsc re- 
gions as "knots" , since they are resolved at 
higher frequcncies. 

• between ~ 10 mas and ~ 30 mas from the 
core, the jet loses its coUimation and be- 

comcs visiblc only in lower rcsolution images 
(Fig. 6, 7). In this region, the jet is oriented 
at a PA of ~ 110° and its opening angle 
is increasing. A well dehned limb bright- 
ened structure is visible (see, e.g., Fig. 7), in 
agreement with the "layercd" structure vis- 
ible in the polarized emission (Aaron 1999). 
As in the hrst region, lower frequcncy images 
reveal a wider jct and a more cvidcnt exter- 
nal shear laycr; this implies that the spectra 
in the inner regions are Hatter than in the 
external regions. 

• at ~30mas from the core the jet shows 

strong bending and a large opening anglc 
(Fig. 8, 9). In this third region the ori- 
entation of the diffuse jet is in agreement 
with the orientation of the kpc-scalc cmis- 
sion. The jet opcning angle is ~43°. The 
radio brightncss is unitorm and there is no 
indication of the hclical structure suggested 
by Conway & Wrobel (1995), which possi- 
bly resulted from the poor (u, -y^-co^erage 
(fivc tclcscopes only) of thcir obscrvations. 
A limb-brightened structurc is visiblc in our 
low resolution images up to ~ 100 mas from 
the core, as can be seen in Fig. 9. 

No radio cmission is dctcctcd cithcr on thc op- 
posite side of the radio core or in the position of 
the putative counterjet feature F reported by Con- 
way & Wrobcl (1995). Howevcr, if wc compare our 
measured total flux densities to single dish data 
(e.g. from the UMRAO database^^), the corre- 
lated VLBI flux density accounts only for a part 

see http : //www . astro . Isa . umich . edu/obs/radiotel/umrao . html 
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of the total flux dcnsity 70% at 5 GHz and 
~ 65% at 15GHz); this is in agreement with the 
presence of a larger scale structure resolved out 
by thc VLBI data, such as the two sidcd cmission 
shown by Ulvestad et al. (1983) and Cassaro et al. 
(1999). 

3.2. Radio Flux Densities and Yariability 

Single dish monitoring of Markarian 501 shows 
that this source is remarkably stable in the radio 
band. Ycnturi ct al. (2001) monitorcd a sample 
of 23 7— and X— ray loud blazars over a 4 year 
period (1996-2000) with the 32 m Medicina radio 
tclcscopc. From thcir data, Markarian 501 is the 
most quiescent source, with an almost flat light 
curve at both 5 and 8.4 GHz. Moreover, Markar- 
ian 501 is part of thc UMRAO databasc, which 
also shows a lairly steady flux density. During the 
period of our observation, the average flux densi- 
tics arc (1.54±0.12) Jy at 4.8 GHz, (1.53±0.13) Jy 
at 8 GHz and (1.29 ± 0.08) Jy at 14.5 GHz. 

We consider the VLBI core flux density us- 
ing our data at 15 GHz (0.6 x 0.9 mas FWHM, 
RA X Dec) , where we have 8 diS^erent epochs (see 
Table 3). Assuming the map peak as the core 
flux dcnsity, wc scc somc dispcrsion in thc val- 
ues, which range between 360 mJy (1997 Apr 25) 
and 508 mJy (1995 Dec 15). This suggests possi- 
blc corc variability, which is dilutcd in singlc dish 
observations, since the VLBI core makes up only 

30% of the total flux density. We note, how- 
evcr, that our obscrvations arc not homogencous 
and were not aimed to monitor flux density vari- 
ability. 

3.3. Proper Motion 

The overall morphology of the source does not 
present dramatic changes between epochs. No 
compact jct componcnt (knot) is prcscnt in the 
jet at all frequencies: if we consider the high- 
est resolution data, the local peaks present in 
thc low rcsolution imagcs bccomc clcarly rcsolvcd. 
Theretore, we limit our study to the innermost 
region, where high-resolution images are avail- 
ablc (thc jct brightncss is too low to bc dctcctcd 
at high frequency beyond lOmas from the core). 
In this region, the jet is always transversally re- 
solvcd: it appcars ccntrally pcakcd at low rcsolu- 
tion, but resolved with a limb brightened structure 



at high rcsolution. This implics that wc could bc 
mixing together different components with diA^er- 
ent properties and moving at different velocities. 
For all thcsc rcasons, wc analyscd imagcs at the 
best angular resolution (typically 0.6x0.9mas, 
RA X Dec) and, in model-Atting, we tried to use 
a numbcr of componcnts that did not producc an 
oversimplified representation of the data. In par- 
ticular, to avoid resolution or frequency effects, 
wc uscd only 15GHz data, plus thc 8.4 GHz ob- 
servation at epoch 1998.48. The results obtained 
from model-fitting of visibilities in Difmap are pre- 
sentcd in Tablc 4 and Fig. 13. Thcsc rcsults arc in 
agreement with those obtained by htting Gaussian 
components to the images (using JMFIT in AIPS). 

The data are in general well fitted by a corc and 
four jet components, in agreement with EP02 (not 
surprisingly, as most data are the same). We la- 
belled components according to thc samc notation 
for ease of comparison: Cl to C4 moving from the 
outermost component inward. However, for the 
reasons statcd above, we adoptcd two Gaussians 
for each of the components C2 and Cl at many 
cpochs. This produced a better fit to the data, in 
agreement with thc rcsolved limb-brightened jet 
structure present in these regions. Nominal posi- 
tion errors are very small. The major uncertainty 
is related to the complcx source structure. We es- 
timate a position error of <0.1 mas for C4, C3, C2 
and ~0.2mas for Cl. 

From our analysis we conclude that no proper 
motion is visible in any of the four componcnts. 
Consider the case of component C2: EP02 mea- 
sured a proper motion of (0.6 ± 0.1)c for this com- 
poncnt. Howcvcr, wc notc that at thc position 
of C2 two faint external pcaks arc visiblc in the 
22GHz image (Fig. 3) on citlicr sidc of the lower 
frcqucncy componcnt. Thcy could cvcn havc a 
diffcrcnt spcctrum (Fig. 10). Thc pcak position 
visible at lower resolution in the 15GHz image is 
thcrclore the blend of the limb-brightencd struc- 
ture visible at 22GHz. In the space VLBI image 
at 5 GHz, the component C2 is at a different posi- 
tion with respect to the 15 GHz images: this could 
be due to different spectra for the different compo- 
nents. We conclude that it is not possible to mea- 
sure a proper motion for such a complex structure 
and that any pattern velocity measurement from 
low rcsolution data is not related to the jet bulk 
speed in this source. 
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3.4. Bulk Motion 

Undcr thc assumption of intrinsically syrnmct- 
ric two-sided jets affected by Doppler boosting 
effects, the observational data allow one to esti- 
matc thc jct bulk vclocity. In thc following, we 
present the results derived from a study of the 
jet/counterjet ratio, core dominance and jet ex- 
pansion. 

3.4- 1- Jet/Counterjet Ratio 

Prom the lack of detection of a counter-jet in 

all thc imagcs (§3.1), wc sct a lowcr limit on the 
jet/counterjet ratio. We use the highest ratio be- 
tween the jet brightness and the counterjet upper 
limit. Wc cstimatc this ratio cxcluding thc inncr- 
most region of the jet; this avoids contamination 
from the unresolved central emission. The best 
imagc for this analysis is thc Spacc VLBI obscr- 
vation obtained in Apr 1998 at 1.6 GHz. From 
this image we can derive a jet/counterjet ratio 
R > 1250 (1 rms) at ^^mas from the core. Given 

the well known rclation R = (^r^4^^) and 

assuming a = 0.5 (see §3.5), this value yields 
/3cos^ > 0.89. In turn, this implies /3 >0.89 and 
6» <27°. 

Our deep, low resolution images allow us to also 
put a constraint at large distance from the core. In 
particular, in the 1.6 GHz image from the 2001.3 
data, at ~ 60 mas from the core we still havc R > 
70; although less extrcmc, this valuc is also vcry 
interesting, as it implies that thc jct rcmains at 
least mildly relativistic (/3 cos d > 0.69) even at 
large distances from the core. 

It is not possible to extend this analysis lurther 
with our data. Howcvcr, wc rcmind thc rcadcr of 
the symmetric morphology of the source on the 
kpc scale (§3.1). This implies that the bulk of the 
jet slows to subrelativistic speed going from the 
parsec to the kpc scales. 

3.4-2. Core Dominance 

Givcn the general relation between total and 
core radio power (Peretti et al. 1984; Giovannini 
et al. 1988, 2001; de Ruiter et al. 1990), we can 
use low frcqucncy unboosted data to compute an 
"expected" value for the core power. A compar- 
ison between the expected and measured 5GHz 
core radio power allows to constrain the jet orien- 



tation and vclocity (Giovannini ct al. 2001). We 
note that at 5 GHz self-absorption effects are small 
(see §4.2); however, to take also into account pos- 
siblc corc variability, wc allow thc obscrvcd core 
flux density to vary by a factor of 2. This is a 
conservative assumption, based on the total flux 
dcnsity monitoring and thc rangc of dcrivcd VLBI 
core flux densities (see § 3.2), and, more generally, 
on statistical flux density variability of AGN in the 
radio. 

To estimate the orientation and jet velocity for 
Mkn 501, we have compared its total unboosted 
flux dcnsity (1.81 Jy) at 408 MHz (Picarra ct al. 
1985), with the nuclear flux density at 5 GHz. 
Since Giovannini et al. (2001) used the arcsecond 
core radio power to derive their correlation, we 
adopted the total VLBI correlated flux density at 
5GHz (~800mJy) as the Mkn 501 core flux den- 
sity. Any possible underestimatc in this assump- 
tion is compensated for by the allowed range (a 
factor of 2) in flux density variability. We found 
that Mkn 501 has to be oriented at 6 < 27° 
with f3 > 0.88. We note that a high jet velocity 
(/3 > 0.95) is allowed only in the range 10° < ^ < 
27°. High velocities arc forbi(l(lcn at smaller an- 
gles; the resulting Doppler lactor would require a 
core much brighter than observed. 

3.4.3. Adiabatic Model 

The functional dependence of the jet intensity 
on the jet velocity and radius for an adiabatically 
expanding jct has bccn discusscd in thc casc of 
relativistic motion by Baum et al. (1997). Since 
the model represents a simpliAcation of the real 
situation, thc rcsults should bc considcrcd with 
caution and compared with other observational re- 
sults. This model has enabled uselul constraints 
to be placcd, in agrccmcnt with othcr data, for 
3C 264 (Baum et al. 1997) and other sources, such 
as NGC315 (Cotton et al. 1999) and 3C449 (Fer- 
etti et al. 1999). 

To apply this model we used our deep VSOP 
observations at 1.6 GHz obtained in March 2001 
and in April 1998. In these observations, the jet 
emission is visible on scales up to ~ 150 mas (see, 
e.g., Fig. 9). First, we derived brightncss proAles 
across the jet: wc uscd thc AIPS task SLICE on 
an image reconstructed with a circular beam of 
7mas diameter, where the longest baselines have 
been omitted to increase the signal to noise ratio. 
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In the innermost region (< 25mas from the core), 
we have 8 shces, taken every half beam-width at 
PA 27.2°. After the main bend, we considered 12 
morc sHccs, cach onc bcam apart, at PA —57.1°. 
We stopped our analysis ~ 100 mas from the core; 
further out, the low jet brightness and the presence 
of extended sub-structures do not allow the data 
to be fitted well. 

Using the AIPS task SLFIT, we htted single 
Gaussians to each prohle. This could bc donc un- 
ambiguously in most cases, although in a fcw shces 
some deviation from a purc singlc Gaussian proAle 
was present. However, the difference between the 
area subtended by the prohle and the fit is always 
smaller than 5% and does not affect the analysis. 
We plot the resulting PWHMs and pcak bright- 
nesses in Fig. 14 and 15; the quantities have been 
deconvolved from the CLEAN beam, according to 
the formula given by Killeen et al. (1986). There 
are clearly two regimes; least-squares fits yield a 
power-law of index — 3.2 in the inner jet (< 30 mas 
from thc core) and —1.3 in thc outer regions. In 
the following, we use these best-ht solutions rather 
then the actual data; this prevents us from spuri- 
ous results brought about by random Auctuations. 

We recall the equations given in Baum et al. 
(1997) for the jet surlace brightness I^: 

I, oc (r,t>,)-(^+2)/3)r-^^*+^^/^^£>2+a (predom- 
inantly parallel magnetic field) 

OC (r^.^;^.)-(55+7)/6)^-(75+5)/6)^2+a (p^g. 

dominantly transverse magnetic iield). 

We assume an injection spectral index 6 = 2, as 
suggested from the average a in the imagc, cxclud- 
ing the self-absorbed core (§3.5). It is not easy to 
decide which magnetic field regime is preferable: 
according to the 8.4 GHz iniagcs of Aaron (1999), 
the Mkn 501 jct is highly polarized at the edges 
with the magnctic hcld parallel to the jet flow, 
while in the jet spine the magnetic held is orthog- 
onal to the jet Aow with a lower degree of order. 
Pollack et al. (2003) report a 5 GHz electric vector 
PA distribution consistent with this result but the 
lower resolution renders the results less conclusive. 
Therefore, since the percentage of ordered field 
and the intrinsic oricntation are not wcU known, 
we consider the two cases of pure parallel and per- 
pendicular magnetic field orientation; in a disor- 
dered magnctic ficld, thc rcal rcsult will bc in bc- 
tween the two extreme situations. Pigure 16 shows 



thc dcrivcd trcnds of /3 for an initial Lorcntz factor 
r = 3.2 and P = 15 (A = 0.95, 0.998) with par- 
allel and perpendicular magnetic fields. In each 
plot, wc draw five lines, corresponding to anglcs 
to the line of sight of 5°, 10°, 25° (i.e., in the 
range of values allowed by the jet sidedness and 
core dominance). 

We note that the jet velocity decreases with 
the core distance, more slowly in the perpendic- 
ular case, and on a shortcr scalc in thc paral- 
lel case. Thc lack of detection of the CJ in all 
our imagcs strongly constrains any possible model 
since it implies a relativistic jet also at largc dis- 
tance from the core (> 50 pc). Models with low 
initial Lorentz factor (P < 5) are dehnitely ruled 
out, regardless of the assumed magnetic ficld ori- 
entation, since they disagree with both the ob- 
served limb-brightened structure (§4.1.1) and the 
jet/counterjet ratio (at least in the case of the 
parallel magnetic field); moreover, they require a 
jet deceleration (§ 4.1.1) between the 7— ray region 
and radio jct rcgion that is too strong. Among 
models with a higher initial speed (L = 15), the 
case of the parallel magnetic field requires an ori- 
entation anglc of 25° , at least at > 70 mas from 
the core, to be in agreement with the observed 
jet/counterjet brightness ratio. We cannot exclude 
models with a narrow angle to the line of sight at 
small distances from the core, where solutions are 
similar. In the case of the perpendicular magnetic 
field, the constraints are less severe since we have 
a fast jet even at large distances from the core. 

3.5. Spectral Index 

In Fig. 10, 11 and 12, wc display thc spcctral 
index distribution at various resolutions between 
1.6 and 4.8 GHz and between 15 and 22GHz. Ob- 
scrvational rcsults arc summarizcd in Tablc 5 for 
high resolution images and described in the follow- 
ing subsections. 

3.5.1. Low Frequency Range: 

At low resolution (3.3 x 1.3 mas, PA = 5°), 
we can study the spectral index distribution at 
distanccs from thc corc largcr than 8 10 mas 
(Fig. 10). The limb-brightened jet structure be- 
tween 8 and 20mas from the core is clearly visi- 
ble. The spectrum is always flat in thc inncr spinc 
(medium value a = 0.08) and steep in the bright 
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shcar rcgions, with valucs in the range 0.9 ^ 2.1. 
In the inner region (< 8 mas from the core), the jet 
is well resolved with a steep spectrum in the more 
extcrnal rcgions (0.7 ^ 1.0) but with a puzzling 
region south of the high brightness region where 
the spectral index is a ~ —0.1. This region is 
clcarly visiblc also in thc lowcr rcsolution imagc 
in Edwards et al. (2000) where it seems connected 
to the inner spine flat spectrum region, while no 
such connection is present in the total intensity 
images. 

Using the space VLBI baselines, we can study 
the sourcc details at higher angular rcsolution 
(1.7 X 0.9 mas, PA= —25°) even at this frequency 
(Fig. 11). The nuclear spectrum is inverted 
(a^;! = -0.45). Exccpt for C4 (confuscd with 
the core and subject to a large uncertainty), the 
jet sub-components show flat to moderately steep 
spectra: 0.1 < a\\g < 0.6. Regions in between the 
components have steeper spectra. 

3.5.2. High Frequency Range: 

At thcsc high frcqucncics, wc arc bcyond thc 
self-absorption turnover and the nuclcar rcgion 
shows a flat, non-inverted spectrum (q: = 0.45), 
generally stccpcning along thc jct. Thc rcsolu- 
tion of 0.6x0.9mas (RA x Dec, Fig. 12) shows 
that the inner regions C4, C3, C2 have in the 
range 0.60 ~ 0.80, while rcgions in bctwccn arc 
much steeper (a ~ 2 and higher). The large com- 
ponent Cl is too faint and resolved to be properly 
imagcd at thcsc high frcqucncics. A marginal sug- 
gestion is present of a spectral steepening moving 
from the inner (spine) to the external jet region 
but at this high angular resolution the jet width 
is not well defined. 

3.5.3. Average Nuclear Spectrum 

It is intcrcsting to considcr thc nuclcar flux dcn- 
sity over the whole range of frequency. We are 
aware that the compact core of Mkn 501 presents 
somc variability (scc §3.2); howcvcr, thc large 
number of observations and frequencies available 
allow us to carry out such a study based on aver- 
age values. 

We plot in Fig. 17 the average core peak 
flux density, from the images at resolution of 
0.6 X 0.9mas (scc Tablc 3); at 1.6GIIz, wc liad to 
consider super-resolved images, in order to avoid a 



large contamination from the jct. Thc spcctrum is 
convex, peaking at ~ 8.4 GHz and steepening to- 
ward the high frequency end. Below the turnover 
frcqucncy, thc spectrum is Aattcr than a = 2.5, 
as would be expected for a single self-absorbed 
component. This suggests that different compo- 
ncnts with diffcrent turnover frequencies may be 
blended together. 

4. DISCUSSION 
4.1. Jet Yelocity 

4.1.1. Yelocity Structure 

Limb brightened jet morphology on parsec 
scales is present in some FR I sourccs but also 
in a few high power FR II sources (Giovannini et 
al. 2002). A possible explanation is thc presence 
of vclocity structurc in thc jct. Whcn vicwcd at 
a relativcly largc anglc to thc linc of sight, this 
structure yields different Dopplcr factors for the 
intcrnal spinc and thc cxtcrnal shcar laycr. In par- 
ticular, an inncr high vclocity spinc could be de- 
boosted while a slower external laycr could be less 
strongly dc-boosted, or evcn boosted. Transversc 
jct vclocity structure has becn found by Laing et 
al. (1999, 2002a,b) in 3C31, where the jet velocity 
in outer regions is ~0.7 the inner spine velocity. 
Chiaberge et al. (2000) invoke such a structure to 
account for obscrvational discrepancies between 
FR I and BL-Lac objccts. 

A possible origin of the velocity structure is the 
interaction of the jet with the surrounding medium 
and its consequent cntrainmcnt; thc mass loading 
slows down the jet in the external regions, while 
the inner spine travels unaffected. However, the 
prcscncc of limb brightcncd jcts in both high and 
low power sources (Giovannini et al. 2002) and the 
evidence in M87 that the jet structure is present 
very near to the core cmission (see the image pre- 
sented by Junor et al. 1999), suggest the possi- 
bility that the velocity structure may originate at 
the basc of thc jct. Mcicr (2003) presents a model 
where the inner high velocity spine is produced 
in the black hole region, while the external shear 
moving at lower velocity is produced in a more 
extended region downstream. 

In Mkn 501, it is interesting to note that the jet 

structurc is visiblc at a distancc < 1 mas from thc 
core, suggesting that it originates very close to the 
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base of thc jct. Howcver, we recall that, at z = 
0.034, Imas corresponds to ~ 0.7pc, i.e., > 10"^ 
Schwarzschild radii (see Rieger & Mannheim 2003, 
and rcfcrcnccs thcrcin for a discussion on thc ccn- 
tral black hole mass in Mkn 501) and, therefore, 
we do not have the angular resolution to see this 
jct inncr rcgion. Nevcrthclcss, we notc that thc 
presence of velocity structure so near to the nu- 
clear source implies too strong a deceleration if it 
is complctcly duc to thc jct intcraction with thc 
ISM. We suggest (in agreement with M87 data of 
Junor et al. 1999) that a transverse velocity gradi- 
ent is intrinsic to thc jct structurc (as discusscd by 
Meier 2003) and that a further velocity decrease 
due to mass loading in the jet £rom the ISM is 
likely to bc prcscnt and possibly dominatcs thc 
jet velocity structure at larger distances from the 
core. 

4-.1.2. Bulk Motion and Jet Orientation 

We have strong indications that in Mkn 501 the 
radio jet emission is oriented at a large angle with 
respect to the line of sight relative to that expected 
for a blazar: 

• From thc jct morphology: thc jct is limb- 
brightened in high resolution images at ^ 
Imas from the core. The different bright- 
ness can be relatcd to a different Doppler 
factor only with an orientation angle 9 ^15°. 
Then, we can reproduce the observed ratio 
bctwccn the brightness of the external shear 
and the inner spine with Tspine = 15, assum- 
ing Tiayer ^ 3 (see Tables 6, 7). With these 
values, thc spinc Dopplcr factor ((^spine) is 
lower than the external layer Doppler factor 
{Siayer)- On the contrary, a jet orientation 
at a smallcr angle and/or a lower spine ve- 
locity imply ^s-pine > ^iayer and therefore a 
centrally peaked structure. 

• Prom the core dominance: the core boosting 

is high but not cxtrcmc; high bulk vclocitics 
require, therefore, an orientation angle 9 in 
the range 10°-27° with respect to the line of 
sight. 

• Prom the fit to the trend of the jet brightness 
and FWHM: assuming an adiabatic model 
and consistency with the jet/counterjet 
brightness ratio, we need a radio jet starting 



with a high vclocity (/3 ~ 0.998), and in the 
case of a parallel magnetic held, has to be 
~25° at >70mas. 

These results seem to be in disagreement with 
constraints derived £rom the high frequency (7- 
ray) emission and variability detcctcd in blazars. 
Studying the rapid burst of TeV photons in Mkn 
421, Salvati et al. (1998) determine that a jet with 
a bulk vclocity T > 10 with opcning anglc and ori- 
entation angle ~ 1/r is required. Similarly, in 
Mkn 501 the expected Doppler factor should be 
>^ 10, as obtained from hts of synchrotron self- 
Compton (SSC) models (Katarzyhski et al. 2001; 
Tavecchio et al. 2001). To reconcile the radio and 
high frcqucncy rcsults wc have to assumc that the 
inner region of the jet (0.001-0.03 pc) is moving 
with a Lorentz factor > 10 and oriented at 6 <5°. 
From this rcgion wc havc high cncrgy radiation 
and negligible radio emission. Over the range from 
0.03 to 50 pc (projected distance) the jet orienta- 
tion has to changc from ~ 5° to 25° (or at least 
~ 15° in thc case of perpendicular magnetic field). 

We propose that this large change in the jet 
direction may be gradual and we derivc thc min- 
imum change requested by observational data at 
different distances from the core. As shown in Ta- 
ble 6, a progressive change in the jet direction with 
respect to the line of sight is possible and in agree- 
ment with the observed pc-scale jet properties. 

At present, there is no obvious explanation for 
this change in the jct dircction in thc inner region. 
It is possiblc that it could bc rclated to helical 
type motion (Yillata et al. 1999) (although helical 
motion is not visiblc in our high rcsolution im- 
ages), or to an instability in the jet due to reasons 
that are thus far unknown. It is also important 
to strcss that small changcs in dircction arc not 
uncommon in FR I radio galaxics and arc mag- 
nified when projection cffccts occur. In any case, 
we rccall that bcyond 30 mas from thc corc thc jct 
position angle becomes stable and closely aligned 
with the kpc-scale structure. 

Wc notc that in thc casc of thc pcrpcndicular 
magnetic field the jet velocity is still very high at 
large distances from the core, in agreement with 
thc lack of dctcction of a visiblc countcrjct in all 
published maps of this source and despite the sym- 
metric structure visible in the VLA low resolution 
image (Ulvestad et al. 1983). However, we recall 
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that 6*^15° is rcquircd to justify thc observed 
Umb-brightened structure from Doppler boosting 
alone, even in the case of perpendicular magnetic 
field. 

4.2. Physical Parameters 

Prom the present data we can derive that: 

• The nuclear region shows a self-absorbed 
spectrum with a peak frequency ~ 
8.4 GHz and a pcak flux dcnsity Sm ~ 
0.55 Jy. PoUowing Marscher (1987), we can 
use these data to estimate the magnetic field 
in the nuclear rcgion (from 0.03 to 0.15 pc): 
B = 3.2x10-5^4^^^-2^(1^^)-!. 

assuming 9 (the core angular size) = 0.2 mas 
(0.15 pc) and (5 = 5 (Table 6), we estimate 
a niagnctic ficld B ~ 0.03 gauss. On the 
otlicr hand, hts of SSC models to a high en- 
crgy state observed in April 1997 yield val- 
ucs of tlic magnctic ficld in thc rangc 0.2- 
0.02 gauss in the inner region (<0.03pc) 
(Katarzynski et al. 2001). This suggests a 
constant or a not vcry largc; dccrcasc of the 
magnetic field from the gamma ray emitting 
region to the beginning of the radio jet. 

• The radio spectra of jet sub-components 
show an evident steepening at high frequen- 
cies, suggesting the effect of radiative losses 
in a high magnetic field. 

• In the region where the jet is well resolved, 
the spectrum is flat in the inner spine con- 
firming thc prcscncc of a higher transport 
eAiciency, lower radiative losses, or reacceler- 
ation mechanisms. The transverse steepen- 
ing from thc spinc to thc shear layer regions 
conhrms the larger losses in these regions, 
probably related to the interaction with the 
external medium. 

• It is not straighttorward to dctcrminc the 
equipartition magnetic field (Pacholczyk 
1970) given the complex structure of the jet. 
Wc can obtain a rough cstimatc by assum- 
ing that electrons and heavy particles carry 
the same amount of energy {K = 1), the 
filling factor = 1, and the frequency range 
is between 10 MHz and 100 GHz. This yields 



~ 0.015 gauss in thc C4 and C3 rcgions and 
0.01 gauss in the C2 and Cl regions. 

• Assuming a viewing angle 9 as in Tables 6 
and 7, the observed angular distances cor- 
respond to the following de-projected linear 
distances: 

- The inner limb-brightened region is present 
at < 4 pc from the core 

- The first jet region extends to ~10mas, 
corresponding to ~ 30 pc 

- The VLBI jet is visible up to 170-270 pc 
(assuming 9 = 25° - 15°, respectively) 

- The kpc-scale structurc is symmetric (i.e. 
not relativistic) at < 20 — 30kpc from the 
core 

- The large scale symmetric structure ex- 
tending to ~ 70" corrcsponds to a lincar sizc 
of 120-200 kpc, in agreement with the ex- 
pected linear size of a FR I radio galaxy with 
the same low frequency total radio power. 

5. CONCLUSIONS 

As a result of the many frequencies and reso- 
lutions availablc, wc havc cstablishcd that a limb- 
brightened structure is indeed present in Markar- 
ian 501, beginning in the very inner jet. This has 
important conscqucnccs: (1) it implies thc prcs- 
ence of velocity structure already at very small 
scale, which jet models will have to take into ac- 
count; (2) it suggcsts, togcthcr with othcr argu- 
ments, that the VLBI jet can not be oriented at 
too small a viewing angle. In particular, the de- 
tcction of a vcry long onc-sided jct and the rcsults 
of an adiabatic expansion fit yield values of the an- 
gle of view in the range 15° < 9 < 25°, depending 
also on the magnetic field orientation. 

The bulk velocity is strongly relativistic at 
the beginning {T ~ 15), decreasing to /3 ~ 0.6 
(r ~1.25) at thc cnd of thc VLBI jct, 100 mas 
from the core (parallel B) or to /? ~ 0.99 (r ~ 7) 
in the case of a perpendicular B. The data con- 
sidered here are consistent with a pattern speed 
of zero. Previous reports, based mostly on lower 
frequency observations, have reported low pat- 
tcrn spccds, howevcr thc conipIcx limb-brightened 
structure revealed at 22 GHz makes measurement 
of component speeds problematic at high frequen- 
cies. 
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Finally, thc study of thc spcctral indiccs has 
revealed a self-absorbed compact core and a com- 
plex spectral index distribution in the jet with 
flat and vcry stccp rcgions. A clcar conncction 
is present between limb-brightening in total inten- 
sity and the spectral index images: the inner spine 
has a Aat spcctrum whilc thc jct morc cxtcrnal rc- 
gions show steep spectra. Our derived magnetic 
field strengths are similar to those obtained from 
data at othcr wavclcngths, yiclding valucs of ordcr 
B ~ 0.03 gauss in the core region. The estimated 
equipartition field in the jet region is in the range 
0.015-0.01 gauss. 
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Table 1 
Obseryation List 



Date 


Frec[uency 


Array 


Observiiig Time 




(GHz) 




(miiij 


1995 Apr T'' 


15 


VLBA 


40 


1995 Apr 12^ 


8.3 


VLBA 


15 


1995 Dec IS'^ 


15 


VLBA 


40 


1996 Apr 23'^ 


8.3 


VLBA 


7.5 


1996 Apr 23 


15 


VLBA 


6 


1996 Jun T 


8.1 


VLBA 


6 


1997 Mar 13^^ 


15 


VLBA 


40 


1997 Apr 25^^ 


15 


VLBA 


90 


1997 Apr 25 


22 


VLBA 


90 


1997 May 26^^ 


15 


VLBA 


50 


1997 May 26 


22 


VLBA 


50 


1997 Aug 4 


1.6 


VLBA-SC+GO+HALCA 


420 


1997 Aug 15<i 


15 


VLBA 


40 


1997 Aug 15 


22 


VLBA 


40 


1998 Apr 7 


4.8 


VLBA+EB+HALCA 


780 


1998 Apr 8 


1.6 


VLBA+RO+GO+HALCA 


600 


1998 Jun 24^ 


8.4 


VLBA+GC+GN+KK+ 


20 






MC+ON+WF 




1999 Jul 19^ 


15 


VLBA-SC 


40 


2001 Mar 5 


1.6 


VLBA+GO+RO+HALCA 


480 


2001 Mar 6 


4.8 


VLBA-HN+HALCA 


540 



NOTE.— EB, Effelsberg (Germany) 100 m, GC, Gilcreek (USA) 26 m GN, Green Bank (USA) 20 m; GO, 

Goldstone (USA) DSN 70 m; HN, Hancock VLBA 25 m; KK, Kokee Park (USA) 20 m; MC, Mcdicina 
(Italy) 32 m; ON, Onsala (Sweden) 20 m; RO, Robledo (Spain) DSN 70 m; SC, St. Croix VLBA 25 m; WF, 
Westford (USA) 18 m. 

"^^LBA 2 cm Survey of Compact Radio Sources (Kellermann et al. 1998) 

''RRPID: Radio Reference Frame Image Database (Fey & Charlot 1997) 

'^ICRE: International Celestial Reference Frame (Ma et al. 1998) 

<iVLBA Project BM082 (Marscher 1999) 
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Table 2 

SuMMARY OF VLBA+HALCA Obseryations 



Date 


Frequency 


Obs. Time 


Tracking 


Other 




(GHz) 


(hr) 


Stations 


Telescopes 


1997 Aug 4 


1.6 


7 


NZ (3'*) 


GO (4^») 


1998 Apr 7 


4.8 


13 


RZ (2.7'*+l.l'^), NZ (3.3^^) 


EB 


1998 Apr 8 


1.6 


10 


RZ [3.6''), NZ (4^*) 


GO, RO 


2001 Mar 5 


1.6 


8 


RZ (l''), NZ (2'') 


RO (4^^), GO (7.5'') 


2001 Mar 6 


4.8 


9 


RZ (0.5''), NZ (l'') 





NOTE. — RZ, Robledo tracking station; NZ, Green Bank tracking station; EB, Effelsberg (100 m); RO, 
Robledo (70 m); GO, Goldstone (70 m). 
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Table 3 
Obseryational Parameters 



Date 


Frcqucncy 


HPBW'^ 


Noisc Lcvcl 


Peak Flux Density 




(GHz) 


(mas) 


(mJy/beam) 


(mJy/beam) 


1995 Apr 7 


15 


0.6 X 0.9 


0.20 


485 






1.2 


0.25 


524 


1995 Apr 12 


8.3 


1.2 


0.40 


524 


1995 Dec 15 


15 


0.6 X 0.9 


0.35 


508 






1.2 


0.30 


546 


1996 Apr 23 


8.3 


1.2 


0.25 


521 




15 


0.6 X 0.9 


0.30 


490 






1.2 


0.35 


520 


1996 Jun 7 


8.1 


1.2 


0.40 


535 


1997 Mar 13 


15 


0.6 X 0.9 


0.15 


381 






1.2 


0.20 


410 


1997 Apr 25 


15 


0.6 X 0.9 


0.10 


360 






1 9 


n 1 

U. lO 


OiyU 




22 


0.3 X 0.5 


0.15 


327 


1997 May 26 


15 


0.6 X 0.9 


0.10 


408 






1.2 


0.10 


438 




22 


0.4 X 0.6 


0.20 


334 


1997 Aug 4 


1.6 


1.5 X 2.9 


2.3 


409 


1997 Aug 15 


15 


0.6 X 0.9 


0.20 


455 






1.2 


0.15 


487 




22 


0.3 X 0.5 


0.25 


350 


1998 Apr 7 


4.8 


0.6 X 0.9 


0.40 


465 






1.2 


0.25 


509 


1998 Apr 8 


1.6 


1.5 X 2.9 


0.12 


421 


1998 Jun 24 


8.6 


0.6 X 0.9 


0.25 


477 






1.2 


0.35 


508 


1999 Jul 19 


15 


0.6 X 0.9 


0.20 


446 






1.2 


0.20 


480 


2001 Mar 5 


1.6 


1.5 X 2.9 


0.50 


316 


2001 Mar 6 


4.8 


1.2 


0.50 


428 



^Column 3: circular beam or RA x Dec in PA 0° 
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Table 4 

ReSULTS 0F THE MODEL - FlTTING 







(mas) 


1995.29 


C4 


0.6 




C3 


2.4 




C2 


4.3 




Cl 


7.6 


1995.96 


C4 


0.7 




C3 


2.6 




C2 


4.0 




Cla 


7.4 




Clb 


7.5 


1996.29 


C4 


0.5 




C3 


2.6 




C2 


4.0 




Cla 


7.7 




Clb 


7.8 


1997.21 


C4 


0.6 




C3 


2.6 




C2a 


3.7 




C2b 


3.9 




Cla 


7.7 




Clb 


7.9 


1997.29 


C4 


0.6 




C3 


2.3 




C2a 


3.7 




C2b 


4.0 




Cl 


7.8 


1997.37 


C4 


0.6 




C3 


2.5 




C2a 


3.8 




C2b 


4.1 




Cla 


7.7 




Clb 


7.9 


1997.62 


C4 


0.7 




C3 


2.4 




C2a 


3.6 




C2b 


4.4 




Cla 


7.3 




Clb 


8.0 


1998.48 


C4 


0.8 




C3 


2.7 




C2 


4.7 




Cl 


7.4 


1999.55 


C4 


0.7 




C3 


2.6 




C2 


4.5 




Cla 


7.8 




Clb 


8.9 
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Table 5 
Spectral Index Results 



Iiiia<2,<' 


l)caiii 


PA 


corc 


C4 




r3 


C3(()C2 


C2 


C2(()C1 


Cl 


"4.8 
«22 


1.7 X 0.9 
0.6 X 0.9 


-25° 
0° 


-0.45 
0.45 


-0.2* 
0.80 


0.6* 
2 


0.1 
0.65 


0.9 


0.6 
0.70 


0.7 
2.4 


0.3 

resolved 



NoTE. — Results of the spectral index imaged in Pigures 11, 12. Typical uncertainties are ±0.05. Data 
marked with an asterisk (*) are estimated after core subtraction and have larger uncertainties. 
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Table 6 

Jet Yelocity Structure - Parallel Magnetic Field 





6 


r 

spine 






^layer 


Notes 


(PC) 


n 












0.0001 - <0.03 


4 


15 


15 


7 


7 


7-ray region 


0.03-0.15 


10 


15 


4 


10 


5 


Radio core 


0.15 - 7 


15 


15 


2 


3 


5 


First jet region 


7-20 


15-20 


15 


2 1 


3 


4-3 


Before of large bending 


20-30 


25 


10-3 


1-2 


2 


2.5 


After the large bending 


50 


25 


1.25 


1.8 


1.1 


1.5 


Final VLBI jet region 



NoTE. — Rcore = projected distance from the core 



Table 7 

Jet Yelocity Structure - Perpendicular Magnetic Field 



(PC) 


e 

n 


r . 

^ spme 


^spine 


^layer 


^layer 


Notes 


0.0001 - <0.03 


4 


15 


15 


7 


7 


7-ray region 


0.03-0.15 


10 


15 


4 


10 


5 


Radio core 


0.15 - 7 


15 


15 


2 


3 


5 


First jet region 


7-20 


15 


15 


2 


3 


4 


Beforo of large bending 


20-30 


15 


10 


2.5 


3 


4 


After the large bending 


50 


15 


10 


2.5 


3 


4 


Final VLBI jet region 



NoTE. — Rcore = projccted distance from the core 
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Fig. 1. — (u, ■u)-coverage of the Space VLBI ob- 
servation at 1.6 GHz (A = 18 cni) on 1998 April 8. 
Available separately as .jpg file figure01.jpg 



Fig. 2. — (u, i))-coveragc of thc Spacc VLBI ob- 
servation at 4.8 GHz (A = 6 cm) on 1998 April 7. 
Available separately as .jpg file figure02.jpg 



Maikarian 501 22GHz August1997 
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Fig. 3.— Isocontour level of Mkn 501 at 22 GHz 
from thc August 1997 epoch. Thc HPBW is 
0.54 X 0.35 mas in PA -11°. The noise level 
is 0.3mJy/beam. Contours are drawn at 1, 
1.5 3, 5, 7, 10, 15, 20, 30, 50, 70, 100, 200, 
and 300 mJy/beam; the greyscale range is 0.3- 
20mJy/beam. The limb brightened structure is 
visible even in the innermost region of the jet. 



lUarkarian 501 15 QHz August 1997 
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Fig. 4. — Isocontour level of Mkn 501 at 15 GHz 
from the August 1997 epoch. The HPBW 
is 0.6 X 0.9 mas (RA x Dec). The noise level 
is 0.2mJy/beam. Contours are drawn at 1, 
1.5 3, 5, 10, 15, 20, 30, 50, 70, 100, 200, 
and 400mJy/beam; the greyscale range is 0.3- 
15mJy/beam. The limb brightening is clearly vis- 
iblc in thc (rcsolyccl) Cl rcgion, but no longcr in 
the inner region because of the lower resolution. 

iMaritarian 501 4.8 GHz April 1998 
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Fig. 5. — Gray scale plus isocontour level of Mkn 
501 from thc 5 GHz VSOP observation in April 
1998. The HPBW is 0.6 x 0.9 mas (RA x Dcc). 
Contours are drawn at 1, 1.5 3, 5, 10, 15, 20, 30, 
50, 70, 100, 200, and 400mJy/beam; the greyscale 
range is 2-20mJy/beam. 
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Markarian 501 4.8 GHz April1998 

5R 1 1 r 




25 20 15 10 5 -5 

Reiative Right Ascension (mas) 

Fig. 6. — Isocontour level of Mkn 501 from the 
5GHz VSOP observation in April 1998. The 
HPBW is 2 X Imas in PA -20°. Contours are 
drawn at -3, 1.5, 2, 3, 5, 10, 20, 30, 40, 50, 100, 
300, and 500mJy/beam. 



Markarian501 1.6 GHz April1998 
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Fig. 7. — Gray scale plus isocontour level of Mkn 
501 from the 1.6 GHz VSOP observation in April 
1998. The HPBW is 3 x 1.5 mas at PA -10°. Con- 
tours are drawn at 1, 1.5, 2, 3, 4, 6, 8, 10, 30, 50, 
100, 200, and 400mJy/beam; the greyscale range 
is 0.80-15 mJy/beam. 




Fig. 8. Gray scalc plus isocontour lcvcl of Mkn 
501 from thc 1.6 GHz VSOP obscrvation in April 
1998. Thc HPBW is 3 x 1.5 mas at PA -10°. Con- 
tours arc drawn at 0.5, 1, 1.5, 2, 3, 5, 7, 10, 30, 
50, 100, and 300 mJy/beam; the greyscale range 
is 0-lOmJy/beam. 
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Markarian 501 1.6 GHz April 1998 
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Fig. 9. — Gray scale plus isocontour level of 
Mkn 501 from the 1.6 GHz VSOP observation in 

April 1998. The image is convoIvcd with circular 
beam of 8.5 mas FWHM. The image noise level is 
0.4mJy/beam. Contours are drawn at 1, 1.5, 2, 3, 

5, 7, 10, 30, 50, 100, 300, and 500mJy/beam; the 
greyscalc range is 0.50-8 mJy/beam. 



Markarlan5D1 LDWfr«s(uene^tt>Ktral lndex (I) AugustlS^ 
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Fig. 10. — Low frequency, low resolution spectral 
index image: af l with 3xl.5mas HPBW (PA 
10°). Contours show total intensity levels at 1.5, 3, 
5, 10, 30, 50, 100, and300mJy/beam(atl.6GHz); 
thc grcyscalc is thc spectral index from —0.3 (in- 
verted) to 0.8 (steep). 



Markarian 501 Low frequency spectral index (11) August1998 
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Fig. 11. — Low frequency, high resolution spec- 

tral indcx imagc: afl with 2x1 mas HPBW (PA 
—20°). Contours show total intensity levcls at 
5, 10, 20, 30, 40, 50, 100, 400 mJy/beam (at 
4.8 GHz); the greyscalc is thc spectral index from 
—0.3 (inverted) to 0.8 (steep). 
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Markarian 501 High frequency spectrai index August 1997 
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Fig. 12. — High frequcncy, high resolution spcc- 
tral index image: all with 0.6 x 0.9 mas HPBW 
(RA X Dec). Contours show total intensity levels 
at 5, 10, 15, 20, 50, 100, and 400mJy/bcam (at 
15GHz); the greyscale is the spectral index from 
(flat) to 1.5 (steep), where S^i/) oc u~". 
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Fig. 13. — Component positions from model- 
Ĕtting. Distances are given in milliarcsecond. AU 
points are at 15GHz, with the exception of the 
1998.48 epoch which is at 8.4 GHz. 
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Fig. 14. — Plots of pcak brightness (top) and 
FWHM (bottom) of Gaussian model hts vs. dis- 
tance from the core. Straight lines represent least- 
squares fits. There is a clear change in behavior 
at about 35 mas from the core. 
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Fig. 15. — Plots of peak brightness vs. FWHM 
of Gaussian model fits. Straight hnes represent 
least-squares hts. The change in behavior visible 
in Fig. 14 corresponds to the change in slope at 
FWHM~ 15mas. The power law slopes are —3.2 
and —1.3. 



Fig. 16. — Results of the adiabatic model fit. The 
low initial velocity with parallel magnetic field is 
at top; middle and bottom plots show high initial 
velocities (/?, = 0.998) for a parallel and perpen- 
dicular magnetic field, respectively. Orientation 
anglcs of 5°, 10°, 15°, 20°, 25° are drawn in oach 
plot, with angle increasing from bottom left to top 
right. For illustration purposes, the vertical axis 
in the bottom plot spans only the range 0.8-1. 
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Fig. 17. — Average core spectrum from data at 
resolution 0.6 x 0.9 (HPBW, RA x Dec). Note the 
turnover at 8.4 GHz and the rapid decrease at 
higher frequency. The low frequency part presents 
a slow growth, suggesting the presence of more 
than one single self-absorbed component. 
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This figure "figure01.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0309285vl 



This figure "figure02.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0309285vl 



